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THE EQUATORIAL AIR GLOW ANEi THE IONOSPHERIC 
GEOMAGNETIC ANOMALY 

ABSTRACT 

OGO-4 observations of OI (6300A) emissions have revealed a global pattern 
in the equatorial airglow hitherto undetected from the ground based observations. 

It is seen that the post-sunset emission rate of OI (6300A) is generally very 
asymmetrical with respect to the geomagnetic equator and shows no apparent 

O 

correlation with the ultraviolet airglow (OI 1304 andl356A) and F-region electron 
density measured simultaneously from the same spacecraft. Both the ultraviolet 
airglow and the ion density measured in the altitude region of 450 km follow simi- 
lar latitudinal variations and exhibit the well known properties of the equatorial 

O 

ionospheric anomaly. The asymmetry in the OI (6300A) emission can be attributed 
to the asymmetry in the height of the F-2 maximum as inferred from the height 
of the maximum emission. From correlative studies of the airglow and the 
ionospheric measurements, the mechanisms for the ultraviolet and the 6300°A 
emission are discussed in terms of the processes involving radiative and dissocia- 
tive recombinations. A relationship between molecular oxygen density and the 
integrated OI (6300A) emission rate is derived and the feasibility of using this 
relationship for estimating 0 2 density is discussed. 
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THE EQUATORIAL AIRGLOW AND THE IONOSPHERIC 
GEOMAGNETIC ANOMALY 


INTRODUCTION 

The existence of the tropical arcs was first noted by Barbier and Glaume 
(1960) from ground based scanning photometers at Tamanrasset (Algeria). Their 

Q 

observations revealed that the 6300A emission was very often in the form of an 
east-west arc usually moving from north to south during the course of the night. 
The extension of these observations to lower latitudes revealed a rather striking 

o 

similarity between the 01 (6300A) emission and F-region equatorial anomaly in 
electron density (Barbier, Weil and Glaume, 1961; Barbier and Glaume, 1962). 

o 

The emission intensity of OI (6300A) and F-region critical frequency followed 
a similar latitudinal pattern; i.e., they both had their minimum values at the 
equator and their maximum values at about 15°-20° on either side of the geo- 
magnetic equator. The observational data which formed the basis of this conclu- 
sion was largely obtained from the African longitudes and was very limited at 

O 

best. Notwithstanding, the success of the Barbier formula relating the 6300A 
zenith intensity to the F-region critical frequency and virtual height led to the 
general belief that the intertropical arcs were of global character and they were 
the manifestations of the same phenomena which gave rise to the geomagnetic 
anomaly. 

These generalizations could not be verified in detail because of the lack of 
observing stations in various longitude zones, particularly in the southern latitudes, 
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but a difficulty became apparent when the observational data from different 
stations in the African and Pacific zones were compared. King (1968) using the 
airglow data obtained from the scanning photometers at Tamanrasset (22.8°N, 
5.5°E) and Agadez (17.0°N and 8.0°E) deduced the rate of the movement of the 
arcs towards the equator and found it to be similar to the rate of the movement 
of the crests of the electron density. VanZandt and Peterson (1969) using the 
observational data from Maui (20.7°N and 156. 3W) obtained a much more complex 
picture of the equatorial airglow than one would infer from the variations of 
critical frequency. They concluded that movement of the F-layer, rather than the 
vertical changes in electron density, is the controlling factor in determining the 
morphology of the equatorial airglow. 

In recent years, the nature of the intertropical arcs has also been studied in 
the ultraviolet region (Hicks and Chubb, 1970; Barth and Schaffner, 1970), par- 
ticularly the 01 lines at 1304 and 1356A, and this has added a new dimension to the 
study of the equatorial airglow. A number of suggestions have been made for the 
possible origin of the ultraviolet airglow in the equatorial region, ranging from the 
chemical processes involving radiative recombination to the direct excitation of 
atomic oxygen by the precipitating energetic particles (Hanson, 1969; Hicks and 
Chubb, 1970; Barth and Schaffner, 1970; Knudsen, 1970; Olson, Peterson and 
Mosley, 1971; Zipf and Stone, 1971). 

The OGO-4 satellite has provided a unique opportunity to study the nature of 
the equatorial airglow since the airglow measurements both in the visible and the 
ultraviolet region were made almost simultaneously with the measurement 
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of a number of ionospheric parameters needed for the understanding of the possible 
physical mechanisms. The purpose of this paper is to present a detailed global 
picture of the equatorial airglow as obtained from the OGO-4 satellite and discuss 
its characteristics in relation to the ionospheric data obtained from the same 
spacecraft. 

Instrumentation 

The OGO-4 satellite was an attitude-stabilized polar orbiting observatory 
covering an altitude range of 400-900 km. The 6300A airglow measurements 
were made using photometers located on the main body and on the OPEP (Orbital 
Plane Experiment Package). The main body photometer was designed to measure 
the integrated vertical emission below the spacecraft. The OPEP photometer, 
through its scan across the horizon, was capable of measuring the vertical 
airglow profiles (Reed & Blamont, 1967, 1972). The ultraviolet photometers 
consisted of gas photoionization chambers sensitive in various regions between 
lOOoA and 1600A. These emissions were identified as 1304A and 1356A lines of 
atomic oxygen as a result of simultaneous observations by. an ultraviolet spep- 
trometer, also on OGO-4. (Hicks and Chubb, 1970; Barth and Schaffner, 1970). 

The retarding potential analyzer aboard the same spacecraft was designed to 
make in-situ measurements of electron and ion temperatures, ion composition, 
and the flux of suprathermal electrons (Chandra, Troy, Donley and Bourdeau, 

1970). 

In late October, 1967, the perigee of the satellite (418 km) was near the 
equator at a local time of about 2130. This is also the altitude region of the F 2 
peak in the equatorial region. This condition, therefore, permitted the most 
meaningful comparison of the equatorial airglow and the plasma density. 
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Discussion of Observational Data 

O 

Figure 1 is a contour map of the integrated emission rate of 6300A airglow 
constructed from the observational data obtained between October 26-30, 1967, 
corresponding approximately to the post-sunset condition. The abscissa and the 
ordinate represent the geodetic longitude and latitude, and the heavy black line 
corresponds to the position of the magnetic dip equator. This local time (2146) 
shown on this map is a nominal local time. All points on this map were obtained 
within 20 minutes of this time, the differences being due to the inclination of the 
orbit, the range in latitude, and the movement of the plane of the orbit. Because 
of the large variation in the column emission rate ranging from less than 25 
rayleighs to values exceedings 800 rayleighs, the map was constructed by choosing 
the contours at logarithmic intervals. Also plotted in this map are the positions 
of the maximum intensity of ultraviolet airglow (heavy dots) as simultaneously 
observed by Hicks and Chubb (1970). For the purpose of discussion we may 

divide the map into three zones centered around -75°, 0, and 120° longitudes. 
These regions are generally referred to as American, African and Asian zones. 
The general features of this map concerning the nature of the equatorial airglow 
may be summarized as follows: 

O 

1. The global distribution of 6300A airglow intensity is very asymmetrical 
with respect to the geomagnetic equator over a wide longitude zone. During the 
post-sunset condition when the northern arc has covered the entire globe parallel 
to the geomagnetic equator, the southern arc is present only in the African zone 
and is almost absent in the American and Asian zones. 

2. There is no indication of this asymmetry in the ultraviolet airglow. The 
positions of the UV intensity maxima (heavy dots) are uniformly distributed in 
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all the longitude zones and are symmetrical with respect to the geomagnetic 
equator. In the American and Asian zones the difference between the latitudinal 
characteristics of the ultraviolet and 6300A airglow is particularly noticeable 
The 6300A arc is absent in the southern latitudes. These observations indicate 
that although similarities can be inferred between the nature of the tropical arcs 
of 6300A airglow and of ultraviolet airglow, their emissions are controlled by 
different processes. 

Figure 2 shows the variations of 6300A emission rate and 0 + density for 

three selected orbits in the Asian, American and African zones. The local time 

for this measurement is about 2130. The altitude of the satellite, which is also 

the altitude at which 0 + is measured, is shown along with the geomagnetic 

latitude. The 0 + variations in this figure show the well known features of the 

% 

geomagnetic anomaly in all the three zones. The density increases continuously 
on either side of the geomagnetic equator, attains maximum values at about ±20° 
latitude and then drops off towards mid-latitudes. General resemblances between 
the latitudinal variations of 6300A airglow and 0 + density can be inferred only 
from the African zone. Both the Asian and American zones show a highly 
asymmetrical variation of 630oA and a clear maximum is indicated only in the 
northern zone. In the American zones, the differences in the behavior of 630oA 
airglow and 0 + are much more evident in the southern latitudes. When we com- 
pare the 0 + variations with the ultraviolet airglow variations (not shown), we find 
much better correlations in all the longitude zones and almost one-to-one cor- 
respondence between the crest of the geomagnetic anomaly and the position of the 
UV arc. We may, therefore, conclude from these observations that the equatorial 
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ionospheric anomaly is much better manifested in the tropical UV arcs than in 
the tropical red arcs. It is interesting to note that ultraviolet airglow observa- 
tions are very recent. The relationship between the red tropical arcs and the 
geomagnetic anomaly has been discussed for almost a decade and was based 

O 

primarily on observational data from the African zone where the 6300A airglow 
follows the variations of 0 + rather remarkably well. 

o 

Both the 6300A and the ultraviolet airglow emission may be caused either 
by chemical recombination involving molecular and atomic oxygen ions or by 
direct excitations of atomic oxygen by precipaitating particles of different energy 

O 

ranges. The 6300A emission may also be caused by thermal excitations of 
atomic oxygen by electrons at high temperatures. The simultaneous measure- 
ment of these parameters from OGO-4 was useful in choosing among these 
various alternatives. Figure 3 is a typical example of the measurement of the 
electron and ion temperatures (T g , T + ) and flux of the suprathermal electrons. 

The plot in this figure corresponds to a pass in the Asian zone where we have 

O 

noted considerable asymmetry in the 6300A airglow. Figure 3 shows that the 
variations of both T e and T + are relatively small in the equatorial region and 
their values are usually less than 1600°K. In this temperature range, any 

O 

contribution to the 6300A emission caused by thermal excitation is almost 
negligible. The latitudinal variations of the suprathermal electrons also do 
not show any pattern which can be associated with the latitudinal variations of 

o 

either 6300A or ultraviolet airglow. In order to explain the observed ultraviolet 
airglow emission by suprathermal electrons, it is necessary to have electron 
fluxes of the order 10 9 /cm 2 -sec in the energy range of 10-20 eV. Our observations^-, 
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therefore, tend to reject the concept that the suprathermal electrons play an 

* 

appreciable role in causing the ultraviolet airglow in the equatorial region. A 
close correlation between the variations of 0 + ions and ultraviolet airglow, on 
the other hand, tends to favor the radiative recombination of electrons and ions 
as a possible mechanism, at least qualitatively (Hanson, 1969). 

, o 

The lack of correlation between O and 6300A in certain longitude zones 

O 

is not at all surprising once we compare the altitude profiles of the 6300A emis- 
sion rate in the various longitude zones. In general, there is an inverse correla- 
tion between the integrated emission rate and the altitude of the maximum emis- 
sion. In Figure 4, we have shown the relative variations of the two parameters 
for the three zones. The integrated emission rate has been obtained from the 
main body photometer, as in Figure 1 and 2, and the altitude of the maximum 
emission rate has been estimated from the scanning photometer. In our repre- 
sentation, we have distinguished the observational data points from the northern 
and southern latitudes by crosses and circles. The main points to be noted from 
Figure 4 are as follows: 

1. The altitude of the maximum emission varies from about 250-350 km in 
both the Asian and American zones. In the African zone this range is mostly 
between 250-300 km. 

2. The altitude of the maximum is higher in the southern latitudes in both 
the Asian and American zones. In the African zones there is no noticeable 
asymmetry. 

3. In both the Asian and American zones, the relationship between integrated 
emission rate, I 6300 (zenith), and the altitude of the maximum, z (I max), is 


approximately 
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log I = (I max) + constant 

H 


( 1 ) 


where H is in the range of 40-45 km. If we identify H with the scale height of 
molecular oxygen, we may infer a neutral temperature T n , in the range of 1300- 
1400°K from the relation H = kT /M(0 2 )gwithk as Boltzmann constant, M(0 2 ), 
the mean molecular mass of 0 2 , and g, the acceleration due to gravity. This 
temperature is characteristic of the high solar flux measured during this period 
and is in reasonable agreement with the neutral temperature inferred from ’the 
measurement of the ion temperature T + (Fig. 3), which is typically in the range 
of 1400-1500°K. A difference of about 100°K between T + and T n is to be ex- 
pected due to the neutral wind heating of the ions (Stubbe and Chandra, 1971). 


Theoretical Discussion 

The relationship between the integrated emission rate and the height of the 
peak emission discussed in the previous section is a direct consequence of the 
dissociative recombination of 0 2 and electrons. Thus the longitudinal asymmetry 
of 630 oA airglow and its lack of correlation with electron and ion density merely 
reflect the zonal asymmetry in the variation of the F 2 peak. The physical 
processes affecting the F 2 peak are (1) plasma diffusion, (2) electric fields, 
and (3) neutral winds associated with the meridional and zonal pressure gradients. 

These processes are usually not independent of each other and have a large 
number of geophysical parameters in common. The theoretical description of 
these processes, therefore, depends upon solving a large number of basic 
ionospheric and atmospheric equations involving time dependent equations of 
continuity, momentum and energy balance. The understanding of the longitudinal 



asymmetry, in addition, would require a three dimensional solution of this problem 
taking into account the variations of the geomagnetic field with respect to the 
geographic coordinates. Notwithstanding the mathematical and numerical com- 
plexity in obtaining such solutions, we may under stand the nature of the basic 
physical processes influencing the characteristics of the equatorial airglow. 

The volume emission rate, Q (6300), from dissociative recombination of O* and an 
electron may be estimated from the following relation (Peterson, VanZant and 
Norton, 1966; Forbes, 1970): 


Q (6300) = 


A (6300) K d y [0 j] N e 
A d + s d (N 2 ) [N 2 ] 


( 2 ) 


where A(6300) and A d are Einstein Coefficients and are .0069 and .0091 sec" 1 , 
respectively. s d (N 2 ) is the rate coefficient for quenching, 7 the rate coefficient 
for the reaction of 0 + ion with 0 2 , and K a the number of excitations per disso- 
ciative recombination. The integrated emission rate I 6300 (zenith) may be 
calculated by integrating Eq. (2) for known distributions of N e , O , and N 2 . 


6300 


(zenith) = 7.6 x 10 -7 K d 7 


r [c 

J 0 1 + 1 .1 x 


[eg N e d z 


10 2 s. (N 2 ) [NJ 


ray 1 ei ghs 


( 3 ) 


Because of the strong quenching, the contribution to I 630Q (zenith) is almost 
negligible below 200 km. Above this altitude, the exospheric temperature becomes 
nearly isothermal and 0 2 and N 2 densities are in diffusive equilibrium and 
may be described by their scale heights. The calculated to I 6300 (zenith) is 
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also very small from the altitude region above 600 km because of the sharp de- 
crease of [0 2 ] . In the altitude region of 200-600 km we may approximate the 
electron density distribution by a simple Chapman function by specifying maxi- 
mum electron density N m , the height of the maximum z m , the scale height of 
atomic oxygen H (O), and a shape factor S. Thus in the altitude region of 200-600 
km, neglecting the height dependence of g, we may write: 


N e (z) = N m expl 


z - z 


1 - 


SH (0) 


- exp - 


(z - z ) 

' m' 

S H (0) 


( 4 ) 


[ 0 ^ = 


to 2 ] z 


exp - 


(z - z ) 

v m ' 

H (0 2 ) 


( 5 ) 


[Nj] ' = [N 2 ] z 

m 


Z - Z ) 

m ' 


h(n 2 


( 6 ) 


where [ 0 2 ] z and Cn 2 ] z are the densities of 0 2 and N 2 at the height of the 

ra m 

F 2 -peak andH(O), H(0 2 ) and H(N 2 ) the scale heights of O, 0 2 and N 2 . From Eqs. 
2-6, we may derive the expressions for I 6300 (zenith) and the height of the maxi- 
mum emission z T (from the condition 3 i/dz = 0)* 

A m ax 

I 630 o (zenith) = 7.6 x 10" 7 K d yN m [0 2 ]; H (0) F (q (zj, S) ( 7 ) 

m 

Z Omax) = 2 m - S H (0) log (1 + 4S - 3.5 s M (z )) (8) 

max 

where F(q (z m ), S) is an integral involving the quenching term, and q (z) and M(z) 
are quenching parameters defined as follows: 



12 


F (q (z m ), S) = S /2l 2 s 



X 2S- 1/2 eX p _ x 

1 + q(* B ) (2 x) 1 - 758 


dx 


1 

x = — exp 
2 


( Z - Z m> 

SH (0) 


q (z) = l.lx 10 2 s d (N 2 ) [N 2 ] z 


M (z) = 


q ( z ) 

l + q ( z ) 


The numerical values of the integral are given for the probable range of q (zj 
and S in Fig. 5 and may be approximated by the following expression: 


F(q(z m ), S) = 


2 . 33S x 2 s r (2 S + 0.5) 


1 + 2 1,75 s q (z m ) 


A r(2S + 0,5) 
s r (.25 S + 0.5) 


B s ~ 

1 + <1 ( Z m) 


with the values of A s and B s given as follows: 


s 

1 

.9 

.8 

.7 

.6 

.5 

.4 

.3 

.2 

.1 

A s 

1.21 

1.20 

1.20 

1.20 

1.20 

1.21 

1.22 

1.24 

1.26 

1.28 

B s 

1.71 

1.41 

1.55 

.93 

.73 

.55 

.40 

.29 

.23 

.29 


From Eq. (8), it is evident that the height of the maximum emission rate is 
always less than the height of the F. -max. Since M(z ) varies between 0 and 1 
from no quenching to complete quenching conditions, the difference z m ~z(I max ) 
is a function of z l and varies in the range of SH (O) log (1 + 4S) and 
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SH(O) log (1 + 0.5S). For atypical value of S= .8 and the neutral atmosphere model 
(Jacchia, 1971) corresponding to T = 1300°K, we estimate z m = 295 and 425 km, 
respectively corresponding to z T = 250 and 350 km. Thus a difference of 
100 km in the height of the peak emission (Fig. 4) near the northern and the 
southern arc corresponds to adifference of 130 km in the heights of the F 2 -peak. 
Because of the geomagnetic control of ionization these differences cause a latitudinal 

O 

displacement of the 6300A arc with respect to the crest in electron density 
(Fig. 3). 

Equation (7) shows the relative importance of the ionospheric and atmospheric 

O 

parameters in determining the 6300A emission rate. Because of the exponential 
relation between 0 2 (z m ) and z m , the changes in the height of the F 2 -peak affect 
the emission rate much more strongly than the changes in N ra and T n . Since T n 
does not vary more than 10-20 percent in the equatorial region, the latitudinal 
variation in I.„„ (zenith) is largely determined by the variations in N and z . 

In the African zone, the F 2 peak as inferred from Fig. 4 does not show any 
significant latitudinal variation. Thus N m is the controlling parameter in de- 
termining the latitudinal characteristic and a close correspondence between the 

O 

crests of the electron density and 6300A arcs is to be expected. In the American 
and Asian zones, the height of the F 2 -peak in the southern latitudes is about 130- 
140 km higher then in the northern latitudes. With the scale height of 0 2 in the 
range of 40-45 km, this difference in the F 2 -peak would cause the northern arc 
to be about 20-30 times brighter than the southern arc and would mask any 
variation associated with N m . It is interesting to note that the asymmetry in the 
F 2 -peak only affects the 6300A. emission rate. The emission rates of 1304A arc 
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relatively insensitive to the changes in the F 2 -peak and are affected only by the 
variation in N 2 (Hanson, 1969). 

The integrated emission rate I 6300 (zenith) may be calculated from Eq. 7 
from the measurements of N and z using a suitable atmosphere model. Unfor- 

m m 

tunately, a detailed comparison between the calculated and observed emission 
rates cannot be made because N m and z m could not be measured directly from 
OGO-4. We may, nevertheless, estimate their values from the measured values 
of N e (z) and z , using Eqs. (3) and (7). In the following, we present a numerical 
example to demonstrate the applicability of Eqs. 7 and 8 for the airglow calcula- 
tions. The calculation is based on the neutral atmosphere model (Jacchia,(1971) 
corresponding to exospheric temperature = 1300°K. The values of y , K d , s d (N 2 ) 
adopted in this calculation are as follows: 

7 = 3.4 x 10 -10 T -1/2 cm 3 - sec -1 (Smith and Fouracre, 1968) 
s d (N 2 ) = 7 x 10" 11 cm 3 -sec -1 (Forbes 1970) 

K d = 0.9 

There is considerable uncertainty in the estimate of K d . Its value has been esti- 
mated in the range of 0.1 and 1 (Forbes, 1970; Brown and Steiger, 1971) from the 
ionospheric measurements. We have adopted a value of 0.9 in close agreement 
with the laboratory measurement (Zipf, 1970). 
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Comparing the observed and calculated values of I 630 o (zenith), we find the 
agreement well within the uncertainty of the various parameters (see table). 

These equations, therefore, not only provide a basis for understanding the nature 
of the equatorial airglow but also provide useful means of estimating 0 2 density 
from the simultaneous measurements of N , z , and I „„ (zenith). 

The most likely cause for the asymmetry in the height of the F 2 -peak is the 
vertical movement associated with the neutral wind [ see Rishbeth, 1972 for 
detailed references] . The neutral wind caused by the global pressure gradient 
tends to drag the ions along the geomagnetic field lines. A meridional wind 
blowing from south to north will tend to move the ionization upward to the south 
of the geomagnetic equator and downward to the north [Abur-Robb and Windle, 
1969]. In this situation the F -peak will be higher south of the equator resulting 
in the asymmetrical distribution in the airglow emission rate. The asymmetry 
in the height of the F 2 -peak may also produce some asymmetry in the crests of 
geomagnetic anomaly inferred from electron density measurements. However 
this asymmetry is relatively much smaller and less noticeable than the asymmetry 
in the 6300A emission rate. The fact that the latitudinal asymmetry in the F 2 - 
peak is confined to the Asian and American zones and absent in the African zone 
cannot be explained in a simple way. Certain asymmetry may be ascribed to the 
relative orientation of the geomagnetic equator with respect to the geographic 
equator and the longitudinal dependence of the magnetic declination angle. But 
the geometrical factor alone does not seem to be sufficient to explain the longi- 
tudinal asymmetry. Any quantitative study of this problem would require a better 
understanding of the global pressure distribution which is the driving force for 


the neutral wind. 


17 


Summary and Concluding Remarks 

In this paper we have presented a global picture of the equatorial airglow 
and have discussed some of the features which were, hitherto, undetected from 
the ground based measurements. Simultaneous measurements of the relevant 
ionospheric parameters helped in understanding the complex nature of the 

O 

6300A emission, particularly its longitudinal asymmetry, in terms of the simple 
physical processes. Our study also provided a basis for estimating the 0 2 density 
in the upper atmosphere from the simultaneous measurements of airglow emission 
rate and readily measurable ionospheric parameters. 

We have not been able to provide an explanation for the longitudinal asym- 
metry in the equatorial wind system. An understanding of this problem would 
not only require a better comprehension of the global pressure gradients but 
a three-dimensional solution of the basic ionospheric and atmospheric equations. 
We hope that the results presented in this paper will stimulate further research 
in this area. 
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FIGURE CAPTIONS 

Figure 1 The global map of the equatorial airglow - showing the contours 

of 6300 A emission rate and the positions of the ultraviolet arc (dots) 
with respect to geomagnetic equator. 

O 

Figure 2 Comparison of the latitudinal characteristics of 6300 A emission rate 
and ion density in the Asian, American and African zones. 

Figure 3 The latitudinal variations of electron and ion temperatures (T e and 
T + ) and the flux of the suprathermal electrons. 

Figure 4 The distribution of the height of the maximum emission rate in the 

O 

Asian, American and African zone and its relation with 6300 A emission 
rate. 

Figure 5 Quenching Integral F(q z , S) as a function of quenching parameter 

m 

q and shape parameter S. 
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